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We propose a model based on a DBI action for the unification of dark matter and dark energy. 
This is supported by the results of the study of its background behavior at early and late times, and 
reinforced by the analysis of the evolution of perturbations. We also perform a Bayesian analysis 
to set observational constraints on the parameters of the model using type la SN, CMB shift and 
BAO data. Finally, to complete the study we investigate its kinematics aspects, such as the effective 
equation of state parameter, acceleration parameter and transition redshift. Particularizing those 
parameters for the best fit one appreciates that an effective phantom is preferred. 
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I. INTRODUCTION 

Unified models of the two main components of the uni- 
verse, dark energy and dark matter, represent an interest- 
ing option for explaining the substantial evidence of the 
current acceleration of the universe. On the one hand no 
observational direct evidence of either of them is avail- 
able, so it might well be the case that they do simply 
not exist and there is a bolder explanation for the effects 
we attribute to them (perhaps extra dimensions are the 
answer). On the other hand, if one believes these two 
entities really fill our universe (in a huge joint proportion 
as compared to the other components) it remains to dis- 
cover what their nature is. In a way, finding out that they 
happen to be manifestations of the same fluid would at 
least simplify the problem in the sense one should have to 
care about the fundamentals of a single fluid. At worst, 
if investigations along these lines were able to refute this 
idea that the two components are just one, we would at 
least be able to face the future in the confidence that 
dark energy and dark matter can be treated separately. 

Following the tradition of trying to find an intercon- 
nection between the world of Particle Physics and Cos- 
mology, it is customary to try and view unified dark en- 
ergy models as scalar field scenarios. One possibility is 
to explore the evolutions contained in a given scalar field 
model, this is actually the approach of this paper. We 
consider a scalar field setup and by fixing some of its 
degrees of freedom we obtain an expansionary cosmol- 
ogy which mimics a dark matter dominated background 
at early times and a dark energy dominated one at the 
late stages of its history. On the other hand, a pop- 
ular procedure to find a scalar fields based motivation 
for a given evolution is to start from a given equation 
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of state and then "reconstruct" the corresponding La- 
grangian (by specifying its kinetic term and potential). 
This widely followed approach has its caveats, however, 
as in general the scalar field model one ends up with is in 
fact a richer scenario, and contains other evolutions that 
the original seed. Leaving aside these remarks, the route 
of scalar fields toward unified dark energy scenarios may 
offer interesting possibilities and our efforts in this paper 
go in this direction. 

It has been put suggested that acceleration is cos- 
mological settings might be the manifestation of non- 
pcrturbative features of of some string theory versions 
This idea has gathered quite a lot of attention as it 
could provide an explanation to early time acceleration, 
that is Jnflation^ (seei, Jl ijl, H, 0, i, i, 0, [HI, E E 
[ii,[il,[ii,|Ii[ii,[liM[I!l,l2jand references therein for 
regular papers and [23| for a review). According to this 
description, the inflation could be a mode accounting for 
the position of a D-brane with three spatial dimensions 
rambling (radially) in a ten-dimensional space-time with 
a warped metric. This interpretation seems to have the 
virtue, among others, that it would allow inflation to pro- 
ceed with much steeper potentials than in the standard 
weakly coupled slow roll inflation model. 

So, a somewhat natural question to ask is this one: 
could a DBI model be responsible for the acceleration we 
observe at present? Moreover, given the similarity of the 
DBI Lagrangian to that of the most popular unified dark 
energy model, the Chaplygin gas, could it also offer a 
satisfactory and perhaps even more suitable alternative 
for the unification of dark matter and dark energy? The 
investigations we report in this paper show that is in- 
deed the case . We construct a purely kinetic DBI model 
for the joint description of the two main components of 
the universe with the bonus that it also displays late- 
time phantom behavior without resorting to true phan- 
tom dark energy at all. 

Any conclusions about the capability of our model to 
represent a solid alternative to other dark energy sce- 
narios must ideally be reached from both the theoretical 
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and the observational perspectives. To that end first we 
carry out the computation and interpretation of the hn- 
ear gauge invariant perturbations of the model. After 
that, we perform a thorough analysis of this novel uni- 
fied dark sector scenario using geometric means: specifi- 
cally we use the SNIa, the BAO, and the CMB shift test. 
These combination of tests allows to take into account the 
early, mid and late time behavior of our model, which is 
expected to have its own features as compared to models 
in which dark energy and dark matter are different com- 
ponents. Our analysis is performed in the Bayesian spirit 
and it allows us to identify the best fit and errors, and to 
complete the information obtained with a computation of 
the evidences on different ranges of the parameters and 
constraints on kinematical quantities of interest. 



II. THE MODEL 

Our scenario is that of a four-dimensional spatially flat 
FRW spacctime filled with a non-canonical scalar field of 
DBI type. Using the customary perfect fluid interpreta- 
tion we set 
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where, in principle, / and V are arbitrary functions. Us- 
age of the symbol 7 was originally motivated by its anal- 
o gy to the Lorentz factor of Special Relativity, given that 
^/((/))(/) is interpreted as the proper velocity of the brane 

Assuming for the above fluid a barotropic equation of 
state of the form p = {T — l)p, we get 



2H 



2l_ 
p 



and the conservation equation reads 
p + SHTp = 0. 



(4) 
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In this paper we explore the case in which both / = /g 
are V = Vq constants. The goal is to obtain a purely 
kinetic model as in other unified dark sector models [24 . 
m, m, 113, ill, so that the field depends solely on the 
scale factor; and as a consequence the same holds for the 
effective pressure and energy density. To that end we 
insert 



7- 



1 



/o 



Vo, 
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into the conservation equation ^ and obtain 

7 + 3/oi?70' = 0. 
Upon replacement of the derivative 7, 

7 = /n7^00- 
and using the Eqs. ([3]), ([7]) we arrive at 



+3i7+ ^ 0, 
4> 7 

which has the following first integral: 



Thus, it is possible to write 
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with c an arbitrary integration constant and ao the value 
of the scale factor today, which we fix as ao = 1 . 

On the one hand we have accomplished our goal of ob- 
taining a DBI model which is purely kinetic throughout 
the evolution. On the other hand, the behavior of the 
model obtained is quite appealing. Using pT|) it can be 
seen that at very late times, i.e. in the regime a ^ ao, 
one has p ~ Vb, whereas at early epochs, i.e. for a <C ao, 
one has p ^ instead. Thus, synthesizing, the solu- 
tion found interpolates between a dust and a de Sitter 
model. Of course, for a positive /qVo it is easy to see 
-ff > 0, so the evolution is indeed expansionary. 

This new model represents an alternative description 
for the unification of dark matter and dark energy, and 
as popular models of this sort, it can be linked to a non- 
canonical scalar field model obtained from the following 
Lagrangian [29j 
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with X = (/)^. Interestingly, the background evolution of 
our model is equivalent to that of a model filled with two 
separately conserved components: a conventional Chap- 
lygin gas (Cg) and a cosmological constant (cc). There 
is, however, an important difference with respect to our 
model, as the equivalence is limited to zeroth order cos- 
mology, i.e. to the background evolution. At the level of 
perturbations, differences are bound to appear, as in our 
model this evolution is realized with a single fluid. It is 
actually possible to derive its equation of state 
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and to draw from it the same conclusion about the 
asymptotic behavior of the model. In contrast, if one 
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would like to stick to a description in terms of two com- 
ponents (Cg+cc) one should demand there is an inter- 
action among them so they act like a single fluid, and 
given the degree of arbitrariness involved in the use of 
interaction terms, we will rather focus on the single fluid 
description provided by the DBI picture. 

Nevertheless, one must keep in mind the possibility 
out of the scope of this paper that the universe is made 
of two fluids, the DBI one plus cosmic dust, which ulti- 
mately should be confronted with the model we present 
here. This possibility has been recently explored from an 
asymptotic behaviour perspective in [30j. 



III. LINEAR PERTURBATIONS 



In the synchronous gauge the line element is given by: 



(14) 



where the comoving coordinates are related to the proper 
time t and position r by dr = dt/a, dx = dr/a, and 
hij is the metric perturbation. The scalar mode of hij 
is described by the two fields h{k,T) and r]{k,T) in the 
Fourier space: 



hij{x, t) = d ke 



ik-X 



kikjh + {kikj - -Sij)6T] 



(15) 

with k = kk. The Einstein equations to linear order ^ in 
k-space, expressed in terms of h and 77, are given by the 
following four equations (s^ : 

er)-]^^h' = A^GaHT°, (16) 
k^ii = ATTGa^{p + p)e, (17) 

h" + 2-h' - 2k^r] = -SnGa'^ST^, (18) 
a 

h" + 6rf + 2-{h' + 67/) - 2^7^ = 
a 

-24TrGa^{p + p)a. (19) 

Here, the quantities 6 and a are defined as (p -|- p)9 = 
ik^ST^, ip + p)a ^ -ik,k, - %/3)Sj and ^ = T; ~ 
(5*r^'/3 denotes the traceless component of the tensor Tj. 
In addition, 9 is the divergence of the fiuid's velocity 
9 = ik^Vj and ' means d/dr. 

Let us consider a fluid moving with a small coordinate 
velocity v'^ = dx'^ /dr, then, can be treated as a pertur- 
bation of the same order as the energy density, pressure 
and metric perturbations. Hence, to linear order in the 
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perturbations, the energy-momentum tensor, with van- 
ishing anisotropic shear perturbation , is given by 



T^ = {p + p}v,^-T^, 
T] = {p + 5p)8]. 



(20) 
(21) 
(22) 



For a fluid with equation of state p = wp, the per- 
turbed part of energy-momentum conservation equations 
T'^.^ = in the k-space leads to the equations 

5' = -(1 + w)(^9+j^- m (^^ - S, (23) 



9' = -nil - 3w)9 - 



Sp/Sp , 2 



1 + w 1 + w 



k'6, (24) 



where 6 = 6p/p and Ti. = a! ja = aH = d. Assuming 
strictly adiabatic contributions to the perturbations, the 
speed of sound for the fluid is 



Sp p 



1 



(25) 



" Sp p j'^ c^/o + ' 

and the time variation of w is 

w' = -3n{l + w){c'i-w). (26) 

Hence, inserting these last two equations in (I^H)) and (|24p . 
they become 

6' = -{l + w)f9 + ^)-3n{cl-w)S, (27) 



9' = -n{l-3cl)9 



1 + w 



-k^6. 



(28) 



Besides, using equations PD]) and ([^ we ar- 

rive at 



h" + Uh' + 37i;2(l + 3c2)J = 0. 



(29) 



At early time, when the overall fluid has w w 0, the 
effective fluid perturbations evolve similar to those of or- 
dinary dust with 9 = 9 = Q, a ^ t^^'^ and from Eqs. 
(|23I29P we obtain 



6" + HS- -n^s = 



(30) 



and S = cit~^ + C2t^^^ , where ci and C2 are arbitrary 
integration constants. In this dust dominated era the 
perturbation grows as (5 ~ a showing an initial unstable 
phase, compatible with the observation that the primor- 
dial universe would have tiny perturbations which seeded 
the formation of structures in the universe. Conclusions 
about the clustering capabilities of other cosmic settings 
with DBI fluids have been studied in [33| . 

At late times, we are interested to find the evolution of 
the linear scalar perturbations for any mode k. To this 
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end we write the second order differential equation for 
the density perturbation S, see [13] 

S" + [1 + 6{cl - w)]HS' + [9{cl - wfH^ 

+3(cf - w')H + 3(c2 - w){H' + H^) + c^k^ 
3, 



-3c4{i + w)ne.{3i) 



Taking into account that in the late time regime the scale 



factor behaves as a cx e 



we can calculate Ti 



Ti. ~ d I a = a(x a. 



(32) 



Considering the expression of 7 given by the Eq. pip one 
obtains the late-time expansion in terms of 1/a for p, p, 
w, cl, w' and TL. In this way, replacing these expansions 
in Eqs. (|28l3ip and keeping only the most significant 
terms one gets 

5" + rSHS' + [9(c2 - wfH^ + 2{cl - w)re + c^fc^] ^ 
= -3c2(l + w)H6i 



(33) 
(34) 



From Eqs. p3|) and (p4|) the evolution of the pertur- 
bation becomes mode dependent with the /Ti.^ term, 
and for low energy modes their solutions can be obtained 
assuming a power law dependence of the perturbations 
with the scale factor, ^ oc a" and oc a**. In this case the 
approximate solutions are given by 



^1 , 



,10 



7]_ 

i5' 



(35) 
(36) 



where 9q, 5i and 82 are integration constants while di 
is a function of 6*0, c and Vq. This shows that the cou- 
pling to in Eq. can be neglected for all scales we 
arc interested on. Also we find that the energy density 
perturbation decreases for large cosmological times for 
modes satisfying the condition k^ jli} ^ 1. For high en- 
ergy modes, k'^ jV? ^ 1, Eq. ([55)1 is hke the equation 
of motion of a dissipative mechanical system. This re- 
semblance emerges using the analogy with the classical 
potential problem 



d 



= -13H<5'2, 



where 



(37) 



(38) 



As for any mode k the potential V has a minimum at 
6 = 0, the function inside the square bracket in Eq. ([57]) 
is a Liapunov function and the perturbation decreases 
asymptotically reaching ^ = in the limit t ^ 00. 



IV. OBSERVATIONAL CONSTRAINTS 

In this section we will set constraints on the param- 
eters of the model from a Bayesian perspective. Our 
analysis will use geometrical tests: the SN type la lumi- 
nositytest, the CMB shift test [H, and the BAO 
test [331. As these two last tests involve early universe 
quantities (the sound horizon at decoupling and dragging 
epochs) , one must consider a slightly more general setup 
and include radiation, which must to be conserved inde- 
pendently from the DBI fluid. This way, the Friedmann 
equation turns out to be 



(39) 

In terms of the fractional energy densities and the red- 
shift we have 



^0 



Mr, ^ ■' 



where 
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foVo - 1 

Pro 



(40) 
(41) 
(42) 



The latter are subject to the normalization condition 



n,. = 1. 



(43) 



In addition, the CMB and BAO tests require that we 
identify a combination of parameters which behaves ef- 
fectively as flm in the high energy regime. In our case 
this mimicry is played by f2c- 

As this paper is a first approach to this model, we 
are setting constraints only on f2c and Qf. In contrast, 
we fix a prior for fif,, taking the WMAP 5-year best fit, 
ilb = 0.0432. Using the tests mentioned before in the 
framework of Bayesian statistics, one should minimize 
the corresponding function in order to obtain flc and 
ilf, see Appendix 1X1 

We have used two different compilations for SNIa data: 
ESSENCE il [111, which combines the first resuhs of 
the survey |40| with the results of Riess et al. de- 
tected by HST[4l| and UNION [H, a vast sample which 
brings together 414 SN from 13 independent datasets: 
recent samples (SLS, ESSENCE), old datasets and dis- 
tant supernovae from HST. In the case of the UNION 
sample, the best values obtained are Qc = 0.256^Q'Qi[Q, 
rif = 0.160t|^:i^J and for the ESSENCE sample fl^ = 



0.257l|]:Jjii^ and % 



0.202 



+0.177 
-0.202 



with the corresponding 



68.30% uncertainties. The lines in the upper sides of the 
plots in Fig[l] represent the locations on the parameter 
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space which correspond to Chaplygin gas cases and the 
points in each of the hnes indicate the case with the low- 
est value. From visual inspection one can infer that 
the Chaplygin ! gas is rejected by our model. In con- 
trast, LCDM (the = locations) is not significantly 
excluded, as for a certain range of rj,„, LCDM cases lie 
in the 68.30% likelihood credible interval. All in all our 
model provides better fits. 

With the aim of compensating for the arbitrariness in 
the choice of priors, we explore different priors on flc and 
fi/. In the case of J7c we have the guidance of all the lit- 
erature of constraints on dark energy constraints which 
more or less suggests preferred regions. To take advan- 
tage of this we explore four priors of different lengths, 
all centered at the value He = 0.25. In contrast, to il- 
lustrate the effect of changing the prior on flf, which is 
a new parameter on which we have no previous clues, 
we divide the physically allowed region fi/ G [0.00,1.0] 
into four equal intervals. From Fig. [2] and Tab. |l] one 
can conclude that among the priors considered, the re- 
gion 0/ e [0.00,0.25], Oc G [0.24,0.26] gives the best 



constraints for the parameters. 



V. MODEL KINEMATICS 

As this is a new model it is worth examining it from 
different perspectives, the kinematic one being a specially 
relevant one. We investigate the redshift dependence of 
the effective equation of state parameter, w{z), and de- 
rived quantities such as the acceleration parameter, g(z), 
or the transition redshift, Zt- In order to obtain the be- 
havior of w{z), we use the expression that relates it with 
the Friedman equation [i^, [4J] 



w{z) 



2dlnH 
3 dz 



(l + z)-l 



Ho 
H 



(44) 



n,xi + zf 



In our model it takes the form 



w{z) 



{fir-ii + z)^ - s^a) ^Jnc^{i + z)^ + - sn/' 



3y^f7e^(l + z)6 + %2 f^^(^i ^ _ + + + Q^^^ {1 + z)^ + n/' 



(45) 



Analyzing the dependence of the equation of state pa- 
rameter with the redshift we can see that the current 
observational data in the context of our model restricts 
it to be smaller than -1, w(z) < —1 with ^\ „ > 
for the current time. More precisely, we obtain w{z ~ 



5+0.052 
-0.081' 



^ „ = 0.081 



+0.081 



=0 

0) 



-0.140 

J 074+0 074 d: 



with the 



0.078' dz \z=0 



0) = -1.052 
UNION sample and w{z 

0.074+°;°^^ with the ESSENCE sample, with the uncer- 
tainties corresponding to the 68.30% interval of confi- 
dence. 

In addition, we study the acceleration parameter, q{z) 
which for instance can be expressed as 



Figs. [3]and[3]depict the shape of the equation of state and 
acceleration parameters with their corresponding 68.30% 
and 95.45% errors. From the analysis of the acceleration 
parameter we gather that there is a strong evidence of the 
transition from a deceleration to an acceleration stage. 
For a better insight on this matter, we have inferred the 
redshift at which the transition happens in different ways. 



ndi + zf 



, ^ 1 
w{z) + -. 



(46) 



Since we have an expression for (/(z), 



9W = 



(i-fz) 4n^(i + z)'V 



i^lix + zf 



f72(l 



(i - ^ni + nj - o, + nr{i + zf + ^VL) + ^i(\^zf^ 



(47) 



we can compute q{zt) = to obtain the transition redshift, zt. For the ESSENCE compilation of SNIa 
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0.20 0.22 0.24 0.26 0.28 0.30 0.32 0.20 0.22 0.24 0.26 0.28 0.30 0.32 



(a) UNION (b) ESSENCE 

Figure 1: Credible intervals from the combination of SN+CMB+BAO observations for two different SN compilation samples. 



(a) UNION 



Prior on fic 

Prior on flf — — — 


flc G [0.24,0.26] 


flc G [0.23, 0.27] 


flc G [0.26,0.28] 


flc G [0.25,0.29] 


flf e [0.00, 0.25] 
flf G [0.25,0.50] 
flf G [0.50,0.75] 
flf G [0.75, 1.00] 


4.534 • 10"''^ 
2.060 • 10"^^ 
4.312 ■ 10"*^ 
5.250 • 10"™ 


2.655 • 10"" 
5.970 • 10"^^ 
1.667 • 10"^'' 
6.557 • 10-^° 


7.797 • lO"'^^ 
1.242 • 10"" 
3.760 • 10"'^* 
8.821 • 10"*^ 


2.258 ■ 10"" 
7.289 • lO"'^^ 
4.681 • 10"'^* 
2.975 • 10~™ 


(b) ESSENCE 


— — Prior on flc 
Prior on flf — - — 


flc G [0.24, 0.26] 


flc G [0.23,0.27] 


flc G [0.26, 0.28] 


flc G [0.25,0.29] 


flf G [0.00,0.25] 
flf G [0.25, 0.50] 
flf G [0.50,0.75] 
flf G [0.75, 1.00] 


2.052 • lO"'*'^ 
1.750 • 10"^^ 
3.435 • 10"^'' 
8.433 • lO"'^^ 


1.194 • 10"*'' 
4.628 • 10"*^ 
6.867 • 10~^° 
4.257 • lO"*^^ 


3.359 • 10"*^ 
9.751 • 10"*^ 
9.742 • 10"*** 
2.501 • 10"^^ 


1.003 ■ 10"*^ 
5.798 • 10"*^ 
1.058 • 10"*'' 
4.115 ■ 10"^^ 



Table I: Bayesian evidences for our unified dark energy DBI model from the combination of SN+CMB+BAO observations for 
two different SN compilation samples. 




(a) UNION (b) ESSENCE 



Figure 2: 3D representation of Bayesian evidences for our unified dark energy DBI model from the combination of 
SN+CMB+BAO observations for two different SN compilation samples. 
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Figure 3: Variation of the equation of state parameter with the redshift for two different SN compilations. 
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Figure 4: Variation of the acceleration parameter with the redshift for two different SN compilations. 



we have Zt = 0.766t°-°47, ^nd for the UNION sample 

^ _ n 77S+0.036 

In [45[ another approach to obtain it was proposed. It 
involves expanding the acceleration parameter, q[z), into 
two terms: 



qiz) = €io+ z 



(48) 



z=0 



Under this definition, we get with the UNION sample 
Zt = OM9t°o°o7l and with the ESSENCE sample zt = 



+0.139 
-0.078- 



These results are in good agreement with 



0.650 

the results obtained in |45l. 1461. 147| . 

Yet another parametrization was considered in [48 

q{z) ^qo + qi 



1 



(49) 



where qo = q{z = 0) is the value of the deceleration 
parameter at present and qi is the parameter that con- 
tains the correction in the distant past {q{z) = qa + qi 
for z ^ 0). With this parametrization we get that 
the value of the transition redshift for ESSENCE is 
Zt = 0.674t°;^^^ and for UNION zt = 0.697^!^:^^^. 

As the value of zt obtained directly by the explicit 
equation (j47p is bigger that the approximate one, we in- 
fer that the approximations arc not good enough for ac- 
counting accurately for the tendency of our DBI fluid to 



inducing a phantom stage. The result obtained by our 
procedure tells us that the acceleration-deceleration tran- 
sition happens before than the other definitions allow to 
estimate. 



VI. CONCLUSIONS 

This paper offers a new alternative to the popular mod- 
els which attempt at a unification of the dark matter 
and dark energy components of the Universe. This new 
model stems from a purely kinetic DBI action, and there- 
fore suggest that these non-canonical actions, which have 
been resorted to as a way to model the early acceleration 
in the Universe, can also serve the same purpose for the 
late time acceleration. At the background level the evolu- 
tion of this model we put forward is like that of a cosmo- 
logical scenario filled with a Chaplygin gas a cosmological 
constant, but this identification ceases when we explore 
further features of the models. To begin with this evolu- 
tion is realized with a single fluid (unlike the mentioned 
Chaplygin gas plus cosmological constant combination) 
and the differences emerge in consequence for quantities 
important at the level of perturbations, specifically the 
difference becomes apparent in the speed of sound of the 
model, so this novel scenario demands an analysis of its 
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own. We can summarize the results as follows. At early 
times the divergence of the velocity perturbation is negli- 
gible, whereas the energy density perturbation is a grow- 
ing one, thus signaling the initial unstable phase required 
for the onset of structures. At late times the velocity and 
energy density perturbations decouple, and the latter be- 
comes negligible as the Universe becomes dominated by 
vacuum energy. 

The observational analysis suggests our model presents 
some attractive features which extend its value beyond 
the theoretical perspective. To begin with current con- 
straints show our model is by far better suited to the 
observations that the most popul ar unified dark sector 
model: the Chaplygin gas [2J, [2^. Our results also in- 
dicate a modest preference of our model as compared to 
the LCDM one. Perhaps the most remarkable outcome 
of this observational analysis is that the best fit corre- 
sponds to a phantom behaviour, i.e. the effective equa- 
tion of state parameter Wcs lies at present below the — 1 
line. It must be remember that this behaviour is achieved 
without actually having to resort to a genuine phantom 
component, so we do not have to be concerned with the 
associated instabilities. 

This study, which has been carried out from differ- 
ent relevant angles and the results achieved, convinces us 
that our model represents a worthy model for the unifi- 
cation of the dark sector, reinforces the theoretical inter- 
ests of DBI models by extending the range of interest to 
the late Universe, and suggests the interest of exploring 
generalizations of this model, probably by relaxing the 
assumption of a purely kinetic Lagrangian, as perhaps 
further degrees of freedom would allow and even better 
suitability to astronomical observations. 



The probability density function p{6\d,Ai) of the pa- 
rameters to have values 6 for the data, d, under the as- 
sumption that the true model is M is provided by Bayes' 
theorem [SOl 



p{e\d,M) 



C{d\e,M)n{e,M) 

J c{d\e,M)TT{e,M)de' 



(A2) 



where p{6\d,A4) and n{6,A4) arc the posterior and 
prior probability density functions (pdf) respectively 
jSOl . [5ll . [H, m, [11] ■ The prior pdf encodes all previ- 
ous knowledge about the parameters before the observa- 
tional data have been collected. It can be regarded as 
a subjective procedure, but its use is compulsory in the 
Bayesian framework, which is the approach used in theo- 
retical frameworks where only one particular realization 
of the measurement is available. 

Parameter estimation in the Bayesian framework 
is based on maximizing the posterior pdf p{6\d,A4), 
whereas in a "strict" frequentist approach one just max- 
imizes C{d\9,M). When one uses flat priors in the 
Bayesian approach then the same conclusions are drawn 
from both approaches and then the difference turns to 
be conceptual only [5l|, [H, . If the measured obscrv- 
ables arc independent form each other and Gaussian dis- 
tributed aroimd their true value, d{6), with a covariance 
matrix, C , given by the experimental errors, maximizing 
C is equivalent to minimizing the chi-square function 



x^e) = (d°''' - d(0)) (d°''-^ - d{e) 

and for uncorrclatcd data C ij — Sijaf, 



(A3) 



Appendix A: STATISTICS AND DATA ANALYSIS 



obs 



die) 



obs 



(A4) 



In the context of a given physical model which depends 
on some parameters, besides fixing the "most likely" val- 
ues of the parameters to yield a series of available ob- 
servational data, one should measure the degree of confi- 
dence in the fact that these data were generated by these 
parameters in an estimated interval. 



1. Parameter estimation 

The likelihood function, C{d\9,A4), is defined as the 
unnormalized probability density of measuring the data 
d = {di, . . . , dn} given the model A4 and taking its pa- 
rameters the values 6 = {6*1, . . . , 6^} jsot . 

Despite our aim to keep the discussion in this section as 
general as possible, when we analyze particular datasets 
we will assume, as usual, that the measurements are nor- 
mally distributed around their true value, so that 



C{d\e,M)(xe-^'^^^/^ 



The second step toward constraining parameters satis- 
factorily is to construct credible intervals [5l| which mea- 
sure the degree of confidence that a certain data was gen- 
erated by parameters belonging to the estimated interval. 

In the Bayesian approach, the credible intervals are 
drawn around the maximum likelihood point, which gives 
the best fit parameters. After obtaining it by the mini- 
mization of the X^(^)j the boundaries of the region con- 
taining 100n% of likelihood are determined as the values 
of the parameters for which has increased by a certain 
quantity 



X 



(A5) 



with 



n = l 



t^'^e'*dt 



Au.n 



(Al) 



t^^'^e-tdt 



' 2' 2 



^'2 



(A6) 
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where ^ (^^, ^^'^ incomplete T function [s^, [57| . 

The la and 2a errors of the parameter 9i are given 
by the 68.30% and 95.45% credible interval contours, re- 
spectively. The upper limit is the maximum value of the 
contour and the lower one the minimum one. 



2. Bayesian evidence 

In Bayes' approach the evidence is an employed tool 
which informs about how well the parameters of the 
model fit the data, after doing an averaging over all the 
parameter values that were theoretically plausible before 
the measurement ever took place [H^ . 

Then the Bayes' evidence is calculated as the average 
likelihood of the model over its prior parameter space, 

E{M) ^ j T^{e,M)C{d\e,M)de, (A7) 

where 'n{0,M.) is the model's prior on the set of param- 
eters normalized to unity (i.e. (J tt{9, A4)d9 = 1.) The 
most common choice is the top hat prior, 7r(0,A^) = l/V 
with V = ria^i (^a.maa: - ^a,mm)- In that case One 
rewrites Bayes evidence as 



£{M) = 



V 



c{e)de. 



(A8) 



One important and unavoidable inconvenient of the 
use of the evidence is its dependence on the prior ranges 
chosen for parameters. In this way we have computed 
the evidences corresponding to different prior ranges from 
comparison in order to find the most suitable one for our 
model, see TablH 

Once we arrived at this point, a remark is required. 
The usual situation in cosmology is that one has more 
than one set of statistically independent observational 
data, {d^^-'}, . . . {d^™-*} in order to constrain the parame- 
ters 6] in that case, one can resort to the joint probability 
density function 



p(6/|d(^^n...nd(™\7W) = 

p{9\d^^\M) X ...xp{e\d^"'\M). 



(A9) 



With the latest rule one can generalize conveniently the 
whole discussion above to the situation with more than 
one dataset. 



Appendix B: ERROR PROPAGATION IN 
DERIVED QUANTITIES 

In our results, the parameters have not symmetric er- 
rors. Then we can not use the standard error propagation 
formula and we have to perform a modification in order 
to account for these non-gaussianities, \E^ . In our case, 
the constraints on the parameters are given in the form, 
where S9i^u and SOi^d are positive quantities. 



The estimated error in a quantity depending on them, 
will be given by an upper limit 



and a lower one 



X^(max(A/,„,-A/,0)' (Bl) 



A// 



^(min(A/,„,-A/,0)' 



(B2) 



where 



A/„ = /(... 0, + A0,„, ,...)- f{e) (B3) 



f{0). (B4) 



This error estimation is based on finite differences, 
however it can be refined if the errors are enough small, 
i.e. A6'i_„ = 59i^u and /S.9i^i = 59.ij. In that case one can 
write 



A/„ ~ 5U 



and 



A// ~ 6fi = 



(B5) 



(B6) 



In Gaussian situations, where A9i^u = A9ij = A9i, 
one gets the standard error propagation formula and 
A/„ = Afi. 



Appendix C: OBSERVATIONAL TESTS 
1. CMB test 

The peaks and troughs of acoustic oscillations are sen- 
sitive to the distance to the decoupling epoch. Therefore 
CMB provides a measure of the ratio of angular diameter 
distances to the decoupling epoch divided by the sound 
horizon size at this time, DA{Zi:)/rs{r^,). Since we have 
assumed a flat universe, instead of Da{z), we can use the 
comoving distance 



De{z) 



dz' 

H{^y 



(CI) 



We can determine the ratio Dc{z^)/r{z^) by the "acoustic 
scale", I A, 
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(C2) 



In this case, we use the fitting function of z^, proposed 
in M 



= 1048 
with 



1 + 0.00124 (nth^) 1 + gi (O^/i 

0.0783 -0-238 



2\92 



91 = 



1 + 39.5 (r^fc/i^) 



0.763 



(C3) 
(C4) 



In order to calculate the drag epoch, we use the formula 
put forward in [63 | 



where 



1291 {nch' 



,0.251 



1 + 0.659 {nji^f^'^^ 



1 + bi (nth^ 



hi = 0.313 [Vtch 
and 



2\ -0.419 



1 + 0.607 {^Ich 



2\0.674' 



(CIO) 



(Cll) 



0.560 



32 



si. 81 ■ 



I + 21.1 {VLbh^Y 
and the comoving sound horizon given by 

, , . da 

rs{z) 



(C5) 



(C6) 



with = 2.469 • 10-5/1-2, c = 2.9979 • 10^ (for Tcmb = 
2.725X) and h = 0.72 [s^. 

CMB also gives a measure of the "shift parameter", 
R{z), which is related to Dc by [6l[ 



R{z,) = Jn,H^D,{z,). 



(C7) 



Constructing a vector containing these quantities, v = 
{Ia, R, z^), and using the maximum likelihood values of 
5-year WMAP [H], v^mb = (302.10, 1.710, 1090.04) one 
can compute the corresponding 

xluB = iv.-vr^^){C-')f;^^{v,-vf^^f (C8) 

where ( ^cmb ^j^^ inverse covariant matrix of the 
data. 

This derivation of WMAP distance priors restricts the 
models to test because it requires that we assume an 
certain cosmological scenario [Sg . 



2. BAO test 

There is a dependence between the peak position of 
the Baryon Acoustic Oscillations (BAO) and the ratio of 
Dv{z) to the sound horizon size at drag epoch, rg (^drag), 
at which the baryons were liberated from photons. Dy {z) 
is a effective distance measure related to the comoving 
distance 



Dv{z) = 



Dl{z) 



H{z) 



1/3 



(C9) 



h2 = 0.238 (r^c/i^ 



, 0.223 



(C12) 



Now, taking into account the Gaussian priors at z = 
0.2 and 0.35 from BAO data appearing in [33|, we calcu- 
late 

XBAo = [vi~Vi ){C )y [vj-vj ) (C13) 



with ^ ^ { 
(0.1980,0.1904). 



Dv(0.2) ' Dv{0.35 



|I} and v^AO = 



3. Type la Supernovae 

The reduced observational data usually reports values 
of the distance modulus 

fith{zi) = 5\ogiQ{dL{z;e)) + fiQ (C14) 

with the dimensionless luminosity distance. 



dL{z;e) = (1 + 2) / 
Jo 



Hodz 



H(z;Ho,e) 



(C15) 



Then the function to minimize takes the form 



2 / n^ i^lthiZj;^lo, 6) - ^obs(zj))'^ 

XsnVA^o, f j - 2^ 2 ■ 



(C16) 



where cr^.j are the measurement variances. But there is 
a nuisance parameter, which makes the computation 
of more intensive as this parameter is marginalized 
over. Often is used an alternative to marginalize it which 
consists in maximi zing the likelihood by minimizing 
with respect to [63|. Then one can rewrite the as 



Xsn(^) = ci - 2c2Mo + C3M0 



(C17) 
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C2 



being The difference with respect to marginahzation over fiQ 

is neghgible in our results (43 |. 

ci^Yl (Mobs(^:j) - SlogipdL izj;9)f ^^^^^ 
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